WIS-ONR-2 
22  June  1953 


Second  Virial  Coefficients  of  Gases  Obeying  a Modified  Buckingham 

(Exp  - Six)  Potential* 
by 

William  E.  Rice^  and  Joseph  O.  Hirschfelder 


Naval  Research  Laboratory,  Department  of  Chemistry,  University  of 

Wisconsin,  Madison,  Wisconsin 


ABSTRACT 


The  second  virial  coefficient  has  been  calculated  for  gases  obeying  a 
modified  Buckingham  potential,  (Exp-Six), 


The  results  are  tabulated  over  the  temperature  range  kT/o  from  0.40  to 
406.,  and  for  seven  values  of  the  parameter  oc  , from  12.  to  15.  The 
treatment  was  entirely  classical,  and  no  corrections  for  quantum  effects 
were  made. 

A separate  paper  will  present  the  evaluation  of  the  parameters  for 
specific  substances  and  applications  of  the  results,  together  with  applications 
of  the  calculated  transport  properties  given  in  the  preceeding  paper. 


* This  work  was  supported  by  Contract  N7onr-28511  with  the  Office 
of  Naval  Research. 

Present  address:  Experiment  Incorporated,  Richmond,  Virginia. 
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L INTRODUCTION 


The  equation  of  state  ox  any  gas  may  be  written  in  virial  form, 

PV  , b(T)  CCT) 

rr  • I + — ♦ — ~ + • • ( 


where  B(T)  is  called  the  second  virial  coefficient.  According  to  classical 
statistical  mechanics,  the  second  virial  coefficient  for  molecules  with 
central  forces  is  directly  related  to  the  energy  of  interaction  between  pairs 
of  molecules,  (p  (r)  , by  the  integral 


fw-v  _ — f s -cp(r)/AT 

3±T  r dr  e 


Up  to  the  present  time,  the  second  virial  coefficient  has  been 
evaluated^  for  Leonard.. Jones  type  potentials,  and,  in  a few  special  cases, 
for  potentials  containing  an  exponential  repulsion  term.  The  Buckingham 
potential  is  more  reasonable  from  a theoretical  standpoint  than  is  a potential 
haying  an  inverse-power  repulsion,  and  is  more  flexible  since  it  contains 
three  parameters  instead  of  two.  Therefore,  we  have  compiled  a complete 
table  of  the  second  virial  coefficient  using  the  Buckingham  potential. 

We  have  taken  the  Buckingham  potential  in  the  following  modified 
form  (which  we  shall  henceforth  refer  to  as  ''Exp-Six")  : 


max  • 


An  extensive  table  of  second  virial  coefficients  for  the  Lennard- Jones 
(12  - 6)  potential  is  given  by  R.  B.  Bird  and  E.  L.  Spots,  University 
of  Wisconsin  CM  - 599  (1950). 
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2 

The  constants  6 , r , and  g characterise  the  molecular  Interaction. 

TXX 

Here  € is  the  maximum  energy  of  attraction;  r is  the  separation  at  the 

m 

energy  minimum;  and  oc  is  a measure  of  the  steepness  of  the  repulsive 
potential.  Fig.  1 shows  the  behavior  of  <p  (r)  . The  position  of  aero  potential 
is  (T  ; its  value  as  a function  of  o'  is  given  in  Table  2.  The  function  defined 
n Eq.  (3)  would  pass  through  a maximum  at  r ^ and  then  approach  - <*> 
as  r -*  0 . The  location  and  height  of  the  maximum  for  each  oc  is  given  in 
Table  2. 


Table  2 

Reduced  separation  <T  * at  which  the  Buckingham  potential 

is  sero,  location  (r  *)  and  the  height  of  the  maximum, 

max 

and  separation  H taken  as  the  lower  limit  of  the  numerical 
integration. 


0 L 

a* 

rr 

max 

<PV*maJ 

ft 

12.0 

12.5 

13.0 

13.5 
U.p 

14.5 

15.0 

0.8761  0051 
0.8798  2937 
0.8831  9813 
0.8862  7341 
0 . 889 1 0396 
0.8917  2672 
0.8941  7037 

0.3024  7469 
0.2730  3685 
0.2469  7188 
0.2238  1785 
0.2031  8787 
0. 1847  5586 
0. 1682  4550 

1705.4129 

3518.3635 

7110.0863 

1 4114.558 

2 7585. 192 
5 3170. 123 

10  1221.60 

0.2784 

0.3186 

0.3442 

0.3620 

0.3844 
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Fortunately  the  height  of  the  maximum  in  bo  large  for  all  practical  cases  that 
very  little  error  is  produced  by  the  anomalous  behavior  of  our  modified 
Buckingham  potential  in  the  neighborhood  of  rmax  The  original  Buckingham 
potential  corresponds  to  the  use  of  Eq.  (3)  for  all  separations.  Instead  of  our 
modification  (Eq.  4)  to  avoid  the  anomalous  behavior  for  small  separations, 
Buckingham  and  Corner^  use  Eq.  (3)  for  separations  larger  than  r^  , 

hat  is  r < r £.  °°  and  the  following  expression  for  separations  less  than  r : 
m m 


(p(r)  * 


_€_[  t «0-r/rJ 
l-L/et[  0 i 


(5) 


This  relation  makes  the  potential  unduly  complicated.  We  have  found  that  the 
effect  of  the  maximum  in  the  potential  manifests  itself  only  for  very  high 
temperatures  and  for  small  values  of  the  parameter  oc 

It  is  convenient  to  introduce  the  following  reduced  quantities; 

r*  = r/rm  . (6) 


<fV> 


<P(r) 

€ 


i r c 
i-Wtc  L * 


e 


,r  <r% 

’ mai,  * 


(7) 


<pV*)*  00, 


o £ r < r 


max . 


(8) 


t ' • JtT/e  . 


(9) 


B(T)/h  b * f TT  N r 3 

® ^ HI 


(10) 


(2) 

R.  A.  Buckingham  and  J . Corner,  Proc.  Roy.  Soc.  A189,  118  (1947). 
According  to  private  correspondence  Buckingham  and  Corner  were  primarily 
concerned  with  interpreting  high  energy  molecular  beam  scattering  which 
is  sensitive  to  the  potential  energy  for  small  separations. 
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In  terms  of  (Exp-Six)  potential,  values  of  B*(  oc  , T*)  are  given  in  Table  1 
for  the  range  0.40  < T*  < 400.  and  12.0  ^ 1 15.0  . 

II.  GENERAL  FORMULAS 

The  integral  for  the  second  virial  coefficient  in  reduced  form  is 


B 


'UT'u-i  C d(fV)  -cpWt-  , 
U>T  T*  J r -TT-  e d 


(ID 


Following  the  procedure  of  Buckingham  and  Corner^  , we  have  divided 
the  range  of  integration  into  two  regions  of  r*  separated  by  the  arbitrarily 
chosen  value  of  the  separation,  M(  ol  ) ; 

p*  «p*\ l *3  d<*V'>  V(0 /t*  . 

B («,T  ) ■ T * J r ~ j7~  e dr 


- T*  { 


, | fJ  d<T(r*)  -9  (r*)/r*  , 

r ■ » ■ e dr1 

T J dr*' 

M(«0 


(12) 


Since  M(  a)  is  taken  to  be  almost  equal  to  (7(a),  the  two  regions  of 
integration  correspond  roughly  to  a separation  between  the  large  separation 
where  the  energy  of  attraction  is  dominan*  from  the  small  separations  where 
the  energy  of  repulsion  prevails.  Different  sorts  of  mathematical  techniques 
are  used,  to  evaluate  the  integrals  in  the  two  regions.  We  shall  employ  the 
symJbol*K(  oc  , T*)  for  the  integral  from  sero  to  M(ct  ) in  Eq.  (12)  , and 
P ( ot  , T*)  for  the  integral  from  M(oc)  to  oo  . Then, 


B*(o^T>-K(o^TVr(*,T*). 


(13) 
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Table  1 (Continued) 
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In  the  range  M(  ot ) < r*  < oo  , j *(r*)  | < | and  T*  varies  from  0 . 40 

to  400.  , the  exponential  function,  exp(-  (p  *(r*)/T*)  , can  be  expanded  in 
powers  of  the  argument  and  the  series  converges  rapidly.  This  expansion  is 
used  in  evaluating  r ( c*  » T*)  and  permits  the  removal  of  the  reduced 
temperature  from  under  the  integral  sign; 


re*,?*) 


04) 


where 


.oo 


J«C«) 


.*3 


M(v) 


d(p*Cr*)  r n. 

dPr~['P’lr')] 


dr' 


(15) 


For  12.0  i a.  < 15.0  , the  expansion  converges  sufficiently  rapidly  for 
the  present  purposes  when  M(  ct)  is  the  largest  integer  multiple  of  0.005 
less  than  <r*(  a)  . This  choice  of  M is  convenient  in  preliminary  hand 
computations  because  it  avoided  interpolation  in  the  numerical  tables  of  the 
exponential  function.  The  coefficients  JQ(  u ) have  alternating  signs  so  that 
Eq.  (14)  can  be  rewritten  in  the  form: 

°»  { \L  j 

rc«.T*)*-E  — — L T*-"'1 

n-6  o!  ' (16) 


It  is  in  the  range  of  integration  of  K(  «,  , T*)  that  our  modification  of  the 
Buckingham  potential  enters  into  the  calculations.  K(  oc  , T*)  may  be 
rewritten  as 


K(<T)*-3 


<p#(r*)/ t* 


dr*  *■  M'l  I-  e 


-9*(*Vv* 


(17) 
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Setting  the  potential  equal  to  °°  for  separations  smaller  than  r 


max 


K(ot 


r*  .MW 


c ?Wt* 


] dr’ 


max 


(18) 


(p*’(M')/T* 


]. 


DU.  NUMERICAL  INTEGRATIONS 

In  evaluating  K(  oc  , T*)  , it  is  certainly  satisfactory  to  ignore  the 

quantity  exp  [-  cp  *(r*)/T*J  in  the  integrand  for  separations  smaller  than 

_ 7 

r ; * , where  this  exponential  function  is  equal  to  4 x 10  for  T*  = 400  . 

For  13.0  < <x  < 15.0  , the  function  actually  is  this  small  for  certain 
separations,  r .*(<*)  • In  actual  practice.,. , the  lower  limit  of  numerical 
integration  of  the  complete  integrand  was  Vt  ( oi ) , which  is  a separation 
smaller  than  r^*(<X)  such  that  (M  - X.  ) is  an  integer  multiple  of  0.0032  . 
This  value  of  )t  was  chosen  so  as  to  facilitate  preliminary  hand  computations. 
For  0C  = 12.0  and  12.  5 , the  complete  integrand  was  used  for  the  entire  range 
between  r^*  and  fct  , so  that  it  and  rrriMW,»  are  the  same.  The  separations 

it  are  given  as  a function  of  oi  in  Table  2.  Having  introduced  it  ( ot)  , 

1 1 ! ' 

in  order  to  shorten  the  range  of  numerical  integration,  we  may  write: 


K(*,T*)  = -it*  + M3[i-e 


(19) 


V 
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The  eight-strip  Newton-Cotes  quadrature  formula  of  the  closed  type 
was  used  to  evaluate  the  integral  appearing  in  K(  « , T*)  in  the  following 
temperature  ranges. 

12.  5 i T*  < 400  . Sixteen  integration  strips  covering  the  range  of  separations 

i 

irom  )t  \ <x.)  to  « ) were  taken  for  T*  = 400.  , 200.  , 100.  , 50.  , 25.  , 
and  12.  5 . The  error  arising  in  the  course  of  this  integration  was  estimated  by 
repeating  the  integration  with  strips  of  half  the  usual  width.  The  errors  for 
several  pairs  of  values  of  o c and  T*  both  within  and  out  of  this  temperature 
range  are  shown  in  Table  3.  The  error  taken  to  be  maximum  for  the  present 
J ‘ temperature  range  is  6 x 10*®  . 

I 

l 

I 


Table  3 


Error  of  the  sixteen-strip  integration  of  K(  ct  , T*)  , the  small- separation 


part  of  the  integral  in  the  second  virial  coefficient,  B*(  ot  , T*}. 


I 

! 


J 


1 . 5625  i T ♦ < 8 . Thirty-two  strips  were  taken  for  T*  = 8.  , 6.25,  4.  , 


3.  125  , 2.  , and  1.  5625  . The  exponential  function  in  the  integrand  was 


neglected  through  the  first  16  strips,  with  no  error  resulting  except  at  T*  -8  , 


j 




T^\^ 

12 

.0 

15 

.0 

B* 

error 

B* 

error 

400. 

0. 1444 

8.  7 x 10-7 

12.  5 

0. 3610 

0.000057 

6.  26 

0. 1688 

0.00013 

0.2681 

0.00023 

3.  125 

-0. 1295 

0.00025 

C . 0050 

0.00045 

1. 5625 

-0.8806 

0.00029 

-0. 6578 

0.00057 

0.78125 

-2.937 

0.00019 

-2. 468 

0.00050 

0.  390625 

0.000016 

0.000022 

0. 40 

-10. 347 

-8.953 
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-where  the  maximum  error  (at  a = 12.0)  was  6 x 10"  out  of  a total 
B*- value  of  0.2215  . 

0 . 80  S T*  S 1 . 50  Thirty-two  strips  were  taken  for  T*  =0.80  to  1.50 
by  steps  of  0. 10.  The  exponential  in  the  integrand  was  neglected  in  the 
first  24  strips,  with  a maximum  resulting  error,  at  T*  = 1. 50  and 
a = 12.0  , of  3.  5 x 10"5  out  of  B*  = -0.  9509  . 

0.40  < T*  < 0.75  Sixteen  strips  were  taken  for  T*  * 0.40  to  0.75  by  steps 
of  0.05  . The  exponential  was  neglected  in  the  first  14  strips,  with  a maximum 
re  salting  error  of  0.00033  for  T*  slightly  above  the  range  (0. 78125)  and 
a = 12. C , where  B*  = -2.937. 

In  the  four  temperature  ranges,  the  errors  in  K(  ot  , T*)  which  are 
quoted  above  are  maximum  values,  which  are  taken  to  be  typical  for  the 
respective  ranges  in  estimating  the  total  error  in  B*(  a , T*). 

For  the  range  of  at  mentioned  above,  and  for  all  n , it  was  found  that 

■7 

no  significant  error  (less  than  one  part  in  10  ) is  introduced  in  evaluating 
the  coefficients  (a.  ),  by  neglecting  the  repulsive  part  of  the  potential  for 
separations  larger  than  A ( cc ) = 14(<x  ) + 2.48  . The  integrand  formed  from 
the  attractive  potential  alone  is  then  integrated  analytically  and  Eq.  (15)  becomes: 

I 


-3-4n 


J (oc)  - (-0 


f ,3  d<TCf*)r  _ . 7n 

J r ~[‘t  <r  >]  dr 


(20) 
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The  eight-strip  Newton-Cotes  quadrature  formula  of  the  closed 

type  was  used  to  evaluate  the  integral  in  J (ct)  . Beginning  at  hl(  oc ) , 

32  strips  of  width  0.005,  then  32  strips  of  width  0.01,  and  finally  40  strips 

of  width  0.05  were  taken.  The  greatest  error  for  any  set  of  eight  stripB 
-7 

was  5.4  x 10  in  the  first  eight  strips  for  oc  - 12.0  and  n = 0 , out  of  a 

total  J (12.0)  - 0.9657716  . This  error  was  taken  to  be  typical  for  J (oc  ) ; 
o o 

the  errors  for  other  values  of  n were  too  small  to  be  apparent  in  the  eight  - 

digit  arithmetic  employed  in  these  calculations.  The  corresponding  error  in 
-7  -1 

B*  is  5.4  x 10  T* 

K(  a , T*)  and  Jn<  oc ) having  been  computed  for  various  arguments, 

it  was  necessary  only  to  perform  a summation  to  obtain  B*(  a , T*)  , which 

was  to  be  tabulated,  r ( oc  , T*)  was  evaluated  up  to  n = 10  . For 

T*  « 0.  390625  and  oc.  = 12.0  the  tenth  term  in  f was  larger  than  for  any 

other  set  of  parameters,  -ing  2.4  x 10”  A rough  extrapolation  indicated 

.5 

the  eleventh  term  to  be  about  5 x 10  For  T*  = 0.  78125  and  a = 12.0  , 

.7 

the  tenth  berm  was  1 . ? x 10  Thus,  the  largest  errors  in  B*  introduced 
by  cutting  off  the  summation  at  n = 10  are  found  in  the  lowest  temperature 
range,  where  they  are  an  order  of  magnitude  smaller  than  the  errors  in 
K(  oc.  , T*)  . The  errors  arising  from  the  numerical  integration  employed 
in  evaluating  J^(  oc)  are  likewise  smaller  than  the  errors  in  the  K-function, 
even  where  the  former  errors  are  at  their  largest. 
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IV.  INTERPOLATIONS 

The  numerical  integrations  described  above  afforded  values  of  B* 

( a , T*)  for  T*  = 0.40  to  0.75  by  steps  of  0.05,  0 . 78125,  0 . 80,  to  1 . 50 
by  steps  of  0.  10,  1.5625,  2..  3.125,  4.,  6.25,  8.,  12.5,  25.,  50.,  100., 
-iOO.,  and  400  . For  other  values  of  T*  , various  interpolation  schemes 
were  used  to  determine  B*  , with  the  aim  of  producing  a more  complete 
table,  in  which  an  uncertainty  of  about  one  part  per  thousand  was  desired. 

The  values  of  B*  for  T*  = 300.  were  found  by  linear  interpolation, 
with  log^T*  being  considered  the  argument  of  the  function  B*  . The 
corresponding  error  for  the  Lennard- Jones  (12  - 6)  potential  was  0.067%  ; 
a slightly  greater  computational  error  would  be  expected  for  the  present 
potential,  because  of  the  greater  curvature  of  B * as  a function  of  T*  in 
this  temperature  range. 

The  values  of  B*  for  T*=  90.,  80.,  70.,  60.,  40.,  30.,  and 
20.  were  interpolated  by  Everett's  method,  with  argument  log^T*  . The 
corresponding  error  in  B*(20)  for  the  (12  *6)  potential  was  0.014%  and  the 
largest  of  the  other  errors  was  0.00  56%  for  T*  = 70.  The  resulting  B*(20) 
values  were  averaged  with  those  found  by  a cur  ve-fitring  method  described 
below. 

The  values  of  B*  for  T*  in  the  range  1. 60  to  20.  were  found  by 


fitting  the  formula, 
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to  several  set*  of  known  points  in  the  range  T*  = 0.  78125  to  200.,  excluding 
the  evenly- spaced  points  between  0.80  and  1 . 50  . For  oc  =12  0,  the 

coefficients  k were  computed  for  several  sets  of  known  points,  as  indicated 

i 

in  Table  4.  The  interpolated  values  obtained  from  various  sets  of  starting 
Doints  are  compared  in  the  table  {*  indicates  that  the  interpolated  point  is  in 
an  interval  adjacent  to  the  middle  point  of  the  range  of  known  points;  **  indicates 
the  B*(20)  computed  by  Everett's  method).  The  discrepancies  for  oc  = 15.0 
were  found  to  be  about  the  same  as  those  for  oc  = 12.0  shown  in  the  table. 
Following  a scrutiny  of  Table  4,  it  was  decided  to  use  the  known  set  T*  = 3.  125 
to  50.  for  interpolating  B*  for  T*  = 7.  to  20 . The  known  set  1.5625  to  12.5 
wap  employed  for  interpolating  in  the  range  3.  2 to  6.0  , and  the  known  set 
0.78125  to  6.  25  for  interpolating  in  the  range  1.60  to  3.  1 . The  B*-value 
obtained  in  this  way  should  be  accurate  to  one  in  the  fourth  decimal  place. 

The  values  of  B*  for  T*  = 0.  85  to  1.  55  by  steps  of  0.  10  were  found 
by  curve-fitting  the  formula, 

a 

B*  = £ k^T*)1  , (22) 

L'-2 

to  the  five  points  in  each  of  the  two  known  ranges:  T*  = 0. 8 to  1. 2 by  steps 
of  0.  1,  and  T*  =1.2,  1.3,  1.4,  1.5,  and  1 . 5625  ; and  interpolating  the 
unknown  points  within  these  ranges.  When  the  formula  was  fitted  to  all 
possible  five-point  ranges  including  T*  = 0.95  and  1.25,  for  the  extreme 
values  of  a , each  of  the  unknown  points  was  an  interior  point  of  four 
different  ranges.  The  discrepancies  amohg  the  four  calculated  values  of 
B*  for  each  of  these  pairs  of  values  of  oc  and  T*  were  at  most  3 x 10-^ 
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At  higher  temperatures.  the  error  of  interpolation  was  less  (b  x 10 
discrepancy  between  two  interpolated  values  of  B*  for  T*  = 1.45)  , 
whereas  for  T*  = 0.85,  die  low-temperature  boundary  of  the  range,  the 
greatest  discrepancy  among  three  interpolated  values  ^*g  1 x 10“*  . 

To  summarize  the  computational  errors  present  in  the  final  values 
of  B*(  , T*)  : For  0.40  < T*  5 0.75,  the  error  ( which  arises  in  K ) 

is  0.000?  ; for  0.80  < T*  1 400.  , the  sum  of  the  interpolation  errors  and 
the  errors  in  K(  cc  , T*)  may  be  as  great  as  0.0002;  except  that,  at  T*  = 300., 
the  larger  error  of  interpolation  may  raise  the  total  uncertainty  to  0.0003  . 
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